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Abstract 

We show that we can also apply the Hirota method to some non- 
integrable equations. For this purpose, we consider the extensions of 
the Kadomtsev-Petviashvili (KP) and the Boussinesq (Bo) equations. 
We present several solutions of these equations. 



1 Introduction 

The Hirota direct method is one of the famous method to construct multi- 
soliton solutions of integrable nonlinear partial differential equations. Hirota 
gave the exact solution of the Korteweg-de Vries (KdV) equation for multiple 
collisions of solitons by using the Hirota direct method in 1971 0. In his 
successive articles, he dealt also with many other nonlinear evolution equa- 
tions such as the modified Korteweg-de Vries (mKdV) 0, sine-Gordon (sG) 
[3!, nonlinear Schrodinger (nlS) JI] and Toda lattice (Tl) [S] equations. 

The first step of this method is to transform the nonlinear partial differ- 
ential or difference equation into a quadratic form in dependent variables. 
The new form of the equation is called 'bilinear form'. In the second step, we 
write the bilinear form the equation as a polynomial of a special differential 
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operator, Hirota D-operator. This polynomial of D-operator is called 'Hirota 
bilinear form'. In fact, some equations may not be written in the Hirota 
bilinear form but perhaps in trilinear or multilinear forms [§]. The last step 
of the method is using the finite perturbation expansion in the Hirota bilin- 
ear form. We analyze the coefficients of the perturbation parameter and its 
powers separately. Here the information we gain makes us to reach the exact 
solution of the equation. 

The equations having Hirota bilinear form possesses automatically one- 
and two-soliton solutions. But when we try to construct the three-soliton 
solutions we come across a very restrictive condition. This condition was used 
as a powerful tool to search the integrability of the equations by Hietarinta 
[7j. Hietarinta also used this condition to produce new integrable equations 
in his articles 0, 0, PH, [II]. 

Most of the works dealt with the Hirota direct method is about the inte- 
grable equations. But in this work, we show that the Hirota direct method 
also can be used to find exact solutions of some non-integrable nonlinear 
partial differential equations. For illustration we consider an extension of the 
Kadomtsev-Petviashvili (KP) equation, 

(u t - Quu x + u xxx ) x + 3u yy + au tt + bu ty + cV 2 u = (1) 

where a, b and c are constants and V 2 w = U X1X1 +U X2X2 + ... + u XmXm for Xj, j = 
1,2, ...,m are independent variables. This extension of the KP equation has 
bilinear and Hirota bilinear form so the Hirota direct method is applicable. 
But when we obtain exact solutions, we shall consider the case c = 0. In this 
case the equation turns out to be 

(u t - Quu x + u xxx ) x + 3u yy + au tt + bu ty = (2) 

and we call it as the extended Kadomtsev-Petviashvili (eKP) equation. This 
equation is integrable if a = b 2 /12. The equation with this condition is 
equivalent to the KP equation. We find exact solution of this equation by 
using the Hirota method for all a, b. Another example to this fact is the 
extension of the Boussinesq (Bo) equation which is 

u u - u xx + 3(u 2 ) xx - u xxxx + auyy + bu ty + cV 2 m = 0. (3) 

Similar to the extension of the KP equation, a, b and c are constants. When 
c 7^ 0, we can also apply the Hirota method to this equation since it has 
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bilinear and Hirota bilinear form. But again, when we consider the exact 
solutions, we shall take c = 0. So we deal with the equation 



u u - Uxx + 3(m 2 ) 



xx 



— U 



xxxx 



+ au vv + bu 



ty 







(4) 



which we call the extended Boussinesq (eBo) equation. The eBo equation is 
integrable if a = b 2 /4. Indeed, under this condition, it is equivalent to the Bo 
equation. The Hirota method gives the exact solutions of the eBo equation 
for all a, b. Before passing to the application of the Hirota direct method, 
let us see how this method works. 

1.1 The Hirota Direct Method 

We review the Hirota direct method in four steps by following Hietarinta's ar- 
ticle [12] closely. Let F[u] = F(u,u x ,u t , ...) be a nonlinear partial differential 
equation. 

Step 1: Bilinearization: We transform F[u] to a quadratic form in the depen- 
dent variables by a bilinearizing transformation u = T[f(x,t, ...),g(x,t, ...)]. 
We call this form the bilinear form of F[u\. Note that for some equations we 
may not find such a transformation. 

Step 2: Transformation to the Hirota bilinear form: 

Definition 1.1. Let S : C n — > C be a space of differ entiable functions. Then 
Hirota D-operator D : S x S — > S is defined as 



where mi, i = 1,2, ... are positive integers and x,t, ... are independent vari- 
ables. 

By using some sort of combination of Hirota D-operator, we try to write 
the bilinear form of F[u] as a polynomial of D-operator, say P{D). Let us 
state some propositions and corollaries on P(D) j!2j . 

Proposition 1.2. Let P(D) act on two differentiable functions f(x,t,...) 
and g(x,t, ...). Then we have 



[D^D?\..]{f.g} 



[(d x -d x ,) m i(d t -d t ,) m *...]f(x,t,...) 

xg(x',t',...)\ 



x'=x,t'=t,... (5) 



P(D){f.g} = P(-D){g.f}. 



(6) 
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Corollary 1.3. Let P(D) act on two differentiable functions f(x,t,...) and 
g — 1, then we have 

P(D){f.l} = P(d)f , P(D){l.f} = P(-d)f. (7) 

Proposition 1.4. Let P(D) act on two exponential functions e Bl and e e ' 2 
where 6i = kiX + wrf + ... + ctj with k iy Wi, ctj are constants for i = 1, 2. 
Then we have 

P(D){e e \e d2 } = P(h - k 2 , a x - a 2 )e e ^. (8) 
For a shorter notation, we use P(pi — p 2 ) instead of P(ki — k 2 , ot\ — a 2 ). 

Corollary 1.5. If we have P(D){a.a} = where a is any nonzero constant 
then we have P(0, 0, ...) = 0. 

Definition 1.6. We say that a nonlinear partial differential equation can be 
written in Hirota bilinear form if it is equivalent to 

m 

EW/V^O, r/ = l,...,r (9) 

a,P=l 

for some m,r and linear operators P^{D). The f %, s are new dependent 
variables. 

Remark 1.7. There is no systematic way to write a nonlinear partial dif- 
ferential equation in Hirota bilinear form. 

Remark 1.8. For some nonlinear partial differential equations we may need 
more than one Hirota bilinear equation. 

Step 3: Application of the Hirota perturbation: We substitute the finite per- 
turbation expansions of the differentiable functions f(x,t, ...) and g(x,t, ...) 
which are 

N N 

f{x,t,...) = f + J2 £m fUx,t,...) , g(x,t,...)=g + ^2e m g m (x,t,...) 

m=l m=l 

(10) 

into the Hirota bilinear form. Here /o, go are constants with the condition 
ifo,9o) 7^ (0)0) to avoid the trivial solution. For the sake of applicability of 
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the method we take the functions f m and g m , m — 1, ...,N as exponential 
functions, e is a constant called the perturbation parameter. For instance 
for N = 2, we take 

/ = fo + e/i + ^ 2 /2 , # = 5o + effi + e 2 g 2 (11) 

where /i = e dl + e 6 * 2 for 6i = hx + aitf + ... + a i: i = 1, 2. We decide what 
the other terms of the functions / and g in the process of the method. 

Step 4: Examination of the coefficients of the perturbation parameter e: 
We make the coefficients of the perturbation parameter e and its powers 
appeared in the Hirota perturbation to vanish. From these coefficients we 
obtain the functions f(x,t,...) and g(x,t,...). Hence by using them in the 
bilinearizing transformation u = T[f(x, t, ...), g(x, t, ...)], we find the exact 
solution of F[u]. 



2 Applications of the Hirota Direct Method 

2.1 The Extended Kadomtsev-Petviashvili (EKP) 
Equation 

The extended Kadomtsev-Pethviashvili (eKP) equation is given by 

(u t - Quu x + u xxx ) x + 3u yy + au tt + bu ty + cV 2 w = (12) 

which is constructed by adding the terms auu, but y and V 2 w = u XlXl + 
u X2X2 + ... + u XmXm multiplied by c to the Kadomtsev-Petviashvili (KP) equa- 
tion where a, b and c are constants and Xj, j = 1,2, ...,m are independent 
variables. Now let us apply the Hirota direct method to the eKP equation. 
Step 1. Bilinearization: We use the bilinearizing transformation 

u(x,t,y) = -2d 2 Jogf (13) 

so the bilinear form of eKP is 

ftxf ftfx fxxxxf ^fxfxxx 3f xx ~\~ 3fyyf ^fy 

rn 

+ aff tt - off + bf ty f - bf t f v + cY,Ux ]X J ~ f%) = 0. (14) 
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Step 2. Transformation to the Hirota bilinear form: The Hirota bilinear 
form of eKP is 

m 

P(D){f.f} = (D t D x + D i x + 3D 2 y + aD^ + bD t D y + cJ2Dl J ){fJ} = 0. (15) 

3=1 

Step 3. Application of the Hirota perturbation: Insert / = 1 + Yl n =i £n fn 
into the equation (|T5|) so we have 

P( J D){/./} = P(D){1.1} + eP{D){f\.l + l.A} 

+ ...+e 2Ar P( J D){/ iV .M = 0. (16) 

Step 4: Examination of the coefficients of the perturbation parameter e: We 
make the coefficients of e m , m = 1,2, ...,N appeared invito vanish. Here 
we shall consider only the case N = 3 and N = A. Note that since eKP is 
not integrable except if a = b 2 /12, we call the solutions obtained using the 
Hirota method as the iV-Hirota solution of eKP. 

2.1.1 N = 3, Three-Hirota Solution of EKP 

Here we apply the Hirota direct method by using the anzats which is used 
to construct three-soliton solutions. We take / = 1 + ef \ + e 2 fi + e 3 fs where 
fi = e 01 + e 02 + e 03 with 9i = hx + u{t + ky + Y^Li r v x j + a i f° r ^ = 1,2,3 
and insert it into (|16jl. The coefficient of e° is identically zero since 

P(D){1.1} = 0. (17) 

By the coefficient of e 1 

P(D){l.h + /i.l} = 2P(d){e 01 + e 02 + e 0:i } = (18) 
we have the relation 

m 

P(pi) = + kf + 3l 2 + auj 2 + buih + c r % = ( 19 ) 

3 

for i = 1,2,3. This relation is called as the dispersion relation. Note that 
when c, the coefficient of V 2 m is not zero, we can apply the Hirota direct 
method. But for simplicity, we take c = in the rest of the calculations. In 
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this case 9i turn out to be 9{ = k { x + U{t + Uy + «j, i = 1, 2, 3. From the 
coefficient of e 2 we get 

(3) 

- P(d)f 2 = - k^Ui - uij) + {h - kj) 4 + 3(l t - I,) 2 

i<j 

+ a(ui - ujj) 2 + b(uji - Uj){li - lj)]e ei+d] (20) 

where (3) indicates the summation of all possible combinations of the three 
elements with i < j. Thus f 2 should be of the form 

/a = A(l, 2)e 9l+02 + A(l, 3)e 6l+e * + A(2, 3)e 92+6s (21) 

to satisfy the equation. We insert f 2 into the equation (|20|) so we get A(i,j) 
as 

AC -\ P (Pi~Pj) 

b(uilj + Ujli) + 2auJiLOj + ki(LOj + 4fc|) + kjfa + 4fcf) - 6k 2 k 2 + 6klj 
b(ujilj + ojjli) + 2auJiUjj + /cj(wj + 4/cj) + kjfa + 4fc|) + 6fcf &j + 6ZjZj 

(22) 

where i,j = 1, 2, 3, % < j. From the coefficient of e 3 we obtain 

P(d){fc} = -[A(l, 2)P(p 3 -p 2 - Pl ) + A(l, 3)P(p 2 - Pl - p 3 ) 

+ A(2, 3)P(p! - p 2 - p 3 )}e e ^ e \ (23) 

Hence f 3 is in the form f 3 = Be 9l+d ' 2+0:i where B is found as 

B _ A(l, 2)P( P3 - gi - p 2 ) + 4(1, 3)P(p 2 - pi - p 3 ) + A(2, 3)P(pi - p 2 - Pa) 

-P(Pl+P2+P3) 

(24) 

The coefficient of £ 4 gives us 
e 29 1+ e 2+ e 3 [Bp{p2 + pg) + 2)A(lj 3)p(p2 _ pg)] 

+ e ei+292+e3 [SP(pi + p 3 ) + 2)A(2, 3)P(pi - p 3 )\ 
+ e ei+e2+2e3 [BP(p! + p 2 ) + 3)A(2, 3)P( Pl - p 2 )] = (25) 

which is satisfied when 

B — A(l, 2)A(1, 3)A(2, 3). (26) 
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To be consistent the two expressions an d (ffi))) should be equivalent i.e. 

A{1, 2)P(p 3 - Pl - p 2 ) + A(l, 3)P(p 2 - Pl - p 3 ) + A(2, 3)P( Pl - p 2 - p 3 ) 



P 



P(Pl + P2 + P3) 

^(1,2)^(1,3)^(2,3). 



(27) 



The above equivalence is satisfied when the condition 

P(Pl - P2)P{Pl - PZ)P{P2 ~ P-i)P{Pl +P2+ P3) 
+P(P1 - p 2 )P(Pl + Pz)P{P2 + P3)P(P3 ~ Pl - P2) 
+P(P1 - P?)P{Pl + P2)P{P2 + P3)P(P2 ~ Pl ~ P3) 
+P(p 2 - P3)P(Pl + P2)P{Pl + Pz)P{pi ~ P2 ~ P3) = 



(28) 



holds. This condition which we call three-Hirota solution condition (3HC) 
can also be written as 



^ P(0"lPl + CF2P2 + 0- 3 P3)P(viPl ~ 02P2)P{P2P2 ~ 0"3P3)P(0"lPl ~ GzVz) = °, 

<T r = ±l 

(29) 

r = 1,2, 3. After some simplifications (3HC) turns out to be 

(12a — b 2 )k\k 2 k 3 2k\w 2 w 3 l 2 l 3 + 2/^2^1^3/1/3 + 2k 2 WiW 2 lil 2 

+ 2k 2 k 3 w 2 w 3 l\ + 2kik 3 wi%v 3 ll + 2k 1 k 2 Wi / w 2 ll + 2k 2 k 3 w\l 2 l 3 

+ 2kik 3 w 2 lil 3 + 2k\k 2 w\l\l 2 — 2k\k 3 w 2 w 3 l\l 2 — 2k\k 3 w\w 2 l 2 l 3 (30) 

- 2kik 2 w 2 w 3 l 1 l 3 - 2k 1 k 2 wiw 3 l 2 l 3 - 2k 2 k 3 WiW 3 l 1 l 2 - 2k 2 k 3 WiW 2 l 1 l 3 

2 ,.2,2 7j2„..2;2 7„2„..2 7 2 7j2„..2;2 ;„2„..2 7 2 7„2„..2;2 



fv^HJ^l^ ^l^2^3 ^2^1^3 ^2^3^1 ^3^2^1 ^3^1^2 



0. 



As we see this condition satisfied when a = 6 2 /12 or for some relations 
on fcj, it;, and U which violate the solitonic property of the solution. The 
coefficients of e 5 and e 6 vanish trivially. Let us focus on the condition ()30|) . 
When the relation a = b 2 /12 holds, the eKP equation is integrable. In fact, 
it is transformable to the KP equation by the transformation 

v! — u, x' — x, y' — y, t' — t + py, 

where p = —b/Q = a/ a/3. If a 7^ b 2 /12, eKP is not integrable. In this case, 
there are other relations which provide f)3(Jj) satisfied. Some of them are; 



8 



Case 1. Any ki = 0, % = 1, 2, 3, the rest are different, 
Case 2. ki = u>i, i = 1, 2, 3, 
Case 3. ki — k, % — 1, 2, 3, 
Case 4. cUj — l^i— 1, 2, 3. 



By using any of these cases, we obtain the exact solutions of eKP. 
2.1.2 iV = 4, Four-Hirota Solution of EKP 

Here we apply the Hirota direct method by using the anzats which is used to 
construct four-soliton solutions. We take / = l+ef\+e 2 f2+£ i f^+s A fi where 
f 1 = e 0i +e e2 +e e 3+e e4 wit j 1 q. = kiX+Uit+ky+aa for i = 1,2, 3, 4 and insert 
it into We will only consider the coefficients of e m , m = 1,2,3,4,5 since 
the others vanish identically. By the coefficient of e 1 

P{D){l.h + h.l] = 2P{d){e 9 ' + e 02 + e 9s + e di } = (31) 

we have the dispersion relation 

P( Pi ) = hui + k\ + 3l 2 + aw 2 + b^ih = (32) 

for i = 1,2, 3, 4. From the coefficient of e 2 we get 

- P(d)h = P(pi - P 2)e 9 ^ + P{pi - ps)e ei+e " + P(p! - p 4 )e 6 ^ 

+ P(P2 ~ P3)e 62+0:i + P(p 2 ~ P4)e° 2+04 + P{p 3 - p 4 )e * +e \ (33) 
Thus j'2 should be of the form 

h = A(l, 2)e dl+d2 + A(l, 3)e Mi + A(l, 4)e ei+e4 

+ A{2, 3)e 62+63 + A(2, 4)e 62+e * + A(3, 4)e e * +6 \ (34) 

where A(i,j), i,j = 1, 2, 3, 4 with i < j are obtained as 

p (Pi ~ Pj) 



Mi, 3) 



P(P* + Pj) 

(3(oJilj + ujjk) + 270;^- + ki{ujj + 4fcJ) + hjfa + Akf) - Qkfk 2 + 61^ 
P(u>ilj + ujli) + 2^uj i uj j + ki(ujj + 4/c|) + kj((Ji + 4Af ) + Qkfk 2 + Qljj 

(35) 
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After some simplifications, the coefficient of e 3 gives 

(4) 

- P(d)f 3 = l A (hj) p (Pm -Pi- Pj) + A(i, m)P(p j - p< - p m ) 

i<j<m 

+ A{j, m)P{ Pl - Pj - p m )]e e ^ +6m (36) 

where (4) indicates the summation of all possible combinations of the four 
elements with % < j < m. Hence fa is of the form 

f 3 = 5(1, 2, 3)e ei+02+93 + B(l, 2, 4)e ei+e2+e4 

+ B(l, 3, 4)e 01+e3+0 i + B(2, 3, A)e e ^ +e \ (37) 

We insert fa into (f3T)J) and obtain 

„/. . \ A(i,j)P(p rn -p i -p j ) + A(i,m)P(p j -p i -p m ) + A(j,m)P(p i -p j -p rn ) 
B(i,j,m) = 

P{Pi + Pj + Pm) 

(38) 

for i,j, m — 1, 2, 3, 4 with i < j < m. From the coefficient of e 4 we have 

PpX/4.1 + f3.f1 + h-f2 + /1./3 + L/4} = 2P(d)h + 2P(D){f 1 .f 3 } + P(D){f 2 .f 2 } 

= 0. 

(39) 

The simplifications gives us that we should have 

B(i, j, m) = A(i, j)A(i, m)A(j, m) (40) 

for i, j, m = 1, 2, 3, 4 with i < j < m. To have consistency the equations (|3~K|) 
and ()4())1 should be equivalent. This yields the condition 

Y P{^iPi+^jPj + ^mPm)P{^iPi-^jPj)P{^jPj-^mP m )P{^iPi-^mPm) = 0. 

(41) 



oy=±l 



for i, j, m, r = 1, 2, 3, 4 with i < j < m, which turns out to be 

(12a — b 2 )k 2 k 2 k 2 m 2k 2 WjW m ljl m + 2k 2 WiW m lil m + 2k 2 a w i Wjlilj 

~\~ 2kjk m zVjZv m l^ ~\~ 2kik m zVj i zv m l j ~\~ 2k*ikjZUiZUjl m ~\~ 2kjk m zu^ljly r i 
+ 2kik m Wjlil m + 2kikjW m lilj 2kik m WjW m lilj 2kik m WiWjljl m (42) 
2kikjWjW m ljl m 2kikjWjW m ljl m 2kjk m WiW m lilj 2kjk m WiWjlil m 



k 2 m 2 I 2 - k 2 m 2 J 2 - k 2 w 2 ] 2 - k 2 vi 2 I 2 - k 2 m 2 J 2 - k 2 m 2 J 2 

hj i W m''j K i W 3 L m ^j aJ i i m K j W mh K m, W j L i li m W i i j 
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where i,j,m = 1,2,3,4, i < j < m. Some of the relations except a = b 2 /12 
which make this condition satisfied are; 



Case 1. Any two of ki = 0, i = 1, 2, 3, 4, the rest are different, 
Case 2. ki = Ui, i = 1, 2, 3, 4, 



Case 3. h = l i} i = 1, 2, 3, 4, 



Case 4. Ui = k, i = 1, 2, 3, 4. 



The equation remaining from the coefficient of e 4 is 

-P{d)U = e 9l+dM9 '[B 123 P(p 4 -p 1 -p 2 -p 3 ) + B 124 P(p 3 - Pl -p 2 -p 4 )} 

+ B lu P(P2-Pl-P3-P4)+B 234 P(Pl-P2-P3-P4)+A(l,2)A(3,A)P(Pl+P2-P3-P4) 

+A(1, 3)A(2, A)P( Pl +p 3 -p 2 - Pi )+A(l, 4) A{2, 3)P( Pl +p 4 -p 2 -p 3 )} = 0. 

(43) 

Thus f 4 = Ce dl+e2+ez+ei where C is obtained as 

C = -[A(l, 2)A(3, A)P( Pl +p 2 -p 3 - p 4 ) + A(l, 3)A(2, A)P( Pl + p 3 - p 2 - p 4 ) 
+A(1, 4)A(2, 3)P( Pl +p 4 - P2 - p 3 ) + B(l, 2, 3)P(p 4 - Pl -p 2 - p 3 ) 
+B(1, 2, A)P{p 3 - Pl - P2 - P4 ) + B{1, 3, A)P{p 2 - Pl - P3 - p A ) 

+5(2, 3, 4)P(p! -p 2 -p 3 - p 4 )} j P{pi +p 2 +p 3 + Pi). 

(44) 

By the coefficient of e 5 we have 

2P{d)f 4 + 2P(£>){/!./ 3 } + P(D){f 2 .f 2 } = (45) 

and when we put the expressions that we have found for f iy i = 1, 2, 3, 4 into 
this equation, it gives 

C = A(1, 2)A(1, 3)A{1, A)A(2, 3)A{2, A)A{3, 4). (46) 

To be consistent the equations (jH| an d should be equal to each other. 
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This yields the condition 
P(Pi ~ P2)P(pi ~ Pz)P{P2 ~ Pz)P{p\ + Pa) 







X P(p 2 + P4)P{P3 + Pa)P{Pa - Pi 


-P2 


-P3j 


+ Pypi 


- P2)P[Pl 


-PA)P{P2 -P4)P{Pl +P3) 










X P(p 2 + P3)P(P3 + Pi)P(P3 - Pi 


-P2 




+ P(Pi 


- P3)P{Pl 


-P4)P{P3 -PA)P{Pl +P2) 










X P(p 2 + P3)P{P2 + Pa)P{P2 ~ Pi 


~P3 


-P4j 


+ P{P2 


~ Ps)P(P2 


- Pa)P{P3 ~ Pa)P{Pi + P2) 










x P(p x + p 3 )P(pi + Pa)P(Pi - P2 






~ P(Pl 


~ P2)P(P3 


- Pa)P{Pi + P3)P{P\ + Pa) 










X P(p 2 + P3)P{P2 + Pa)P{Pi + P2 


-P3 


-P4) 


~P(Pl 


- Ps)P(P2 


- p 4 )P(pi + P2)P{P\ + Pa) 










X P(p 2 + P3)P{P3 + Pa)P(Pi + P3 


-P2 


-P4) 


~P{P1 


- Pa)P(P2 


- P3)P{pi + P2)P{pi + Pa) 










x P(p 2 + p 4 )P(p 3 + Pa)P{Pi + Pa 


-P2 


-Pa) 


~ P(Pl 


~ P2)P{P1 


- P3)P{P\ - Pa)P(P2 - Ps) 







x P(p 2 - p 4 )P(p 3 - Pa)P(pi +P2+P3+ Pa) = 0. 

(47) 

which can also be written as 

4 

53 P(5>iF*) I] [P^ iPi -a jPj )} = 0. (48) 

CTi=±l 1=1 0<i<jr'<4 

We call this condition as four-Hirota solution condition (AHC). Here the 
question is whether the cases satisfying (|4*2|l also satisfy {AHC) automati- 
cally. In the hand we know a case which satisfies both conditions which is 

Case 1. Any two of ki = 0, % = 1, 2, 3, 4, the rest are different. 

Now, for illustration, let us see the graphs of the two- and four-Hirota solu- 
tions of eKP. Here in order to get the solutions we give arbitrary values to a, 
b, ki and Wi and from the dispersion relation we obtain U. Note that in our 
choice a ^ b 2 /12. 
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i) N = 2, The Two-Hirota Solution of EKP: 

The constants are 



a = 2, b = 9, k\ = 1, k 2 = 3, 




Figure 3: t=-2 Figure 4: t=2 
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Figure 5: t=4 Figure 6: t=6 

ii) iV = 4, The Four-Hirota Solution of EKP: 

The constants are chosen according to the Case 1 and the dispersion rela- 
tion. The constants are, 

a = 2, b = 9, h = 0, k 2 = 0, k 3 = 1, k 4 = 2, 



wi = 4, w 2 = -2, w 3 = 3, w 4 = -5, 




Figure 7: t=-6 Figure 8: t=-4 
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Figure 9: t=-2 



Figure 10: t=2 




2.2 The Extended Boussinesq (EBo) Equation 

The extended Boussinesq (eBo) equation is given by 

u tt ~ u xx + 3(u 2 ) xx - u xxxx + au y y + bu ty + cV 2 u = (49) 

which is constructed by adding the terms au yy , bu ty and V 2 w = u xlXl + 
Ux 2 x 2 + ••• + Ux m x m multiplied by c to the Boussinesq (Bo) equation where 
a, b and c are constants and Xj, j = 1,2, ...,m are independent variables. 
Now let us apply the Hirota direct method to the eBo equation. Step 1. 
Bilinearization: We use the bilinearizing transformation 

u(x,t,y) = -2d 2 Jogf (50) 
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so the bilinear form of eBo is 



f ftt ft fxxf fx fxxxxf ^fxfxxx 3f xx 

III 

+ afyyf - afl + bft y f - bftf y + cJ2(f^f - 4) = 0. (51) 

i=i 

Step 2. Transformation to the Hirota bilinear form: The Hirota bilinear 
form of eBo is 

rn 

P(D){f.f} = (A 2 - Dl - D x + aD\ + bD t D y + c£ ££,){/ ./} = 0. (52) 

i=i 

Step 3. Application of the Hirota perturbation: Insert / = l+$^ n =i ^"/ninto 
the equation (jo^j) so we have 

P(D){f.f} = P(D){1.1} + eP(D){A.l + l./i} 

+ ...+£ 2Ar P( J D){^.M = 0. (53) 

Step 4: Examination of the coefficients of the perturbation parameter e: We 
make the coefficients of e m , m = 1, 2, N appeared in (|53|) to vanish. Here 
we shall consider only the case N = 3 and N = A. 

2.2.1 N = 3, Three-Hirota Solution of EBo 

Here we apply the Hirota direct method by using the ansatz which is used 
to construct three-soliton solutions. We take / = 1 + ef \ + s 2 f2 + e 3 fs where 
fi = e dl +e 92 + e dz with 6i = kiX + U{t + ky + X^7=i r u x i + a « f° r * = 1' 2, 3 and 
insert it into (|53p. The coefficient of e° is identically zero. By the coefficient 
of e 1 

P(D){l.f x + f x .l} = 2P(d){e dl + e° 2 + e 03 } = (54) 
we have the dispersion relation 

m 

P{ Pi ) = - kf - kf + all + ivik + c Y. r i = ( 55 ) 

3 

for i — 1, 2, 3. Similar to the eKP equation, we see that when c, the coefficient 
of V 2 u is not zero, we can apply the Hirota method. But for simplicity we 
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take c = in the rest of the calculations. In this case 9i become 6{ = 
kiX + u>it + lit) + CKj. From the coefficient of e 2 we get 

(3) 

- P{d)f 2 = - cuj) 2 - (hi - kj) 2 - (ki - k f + a(h - l 3 ) 2 

i<j 

+ b(u i -uj j )(l i -l j )]e ei+(> i (56) 

where (3) indicates the summation of all possible combinations of the three 
elements with i < j. Thus f 2 should be of the form 

h = A(l, 2)e 9l+9 * + A(l, 3)e 9l+es + A(2, 3)e e *+ 03 (57) 

to satisfy the equation. We insert f 2 into the equation (J56j) so we get A(i,j) 
as 

2uiUj — 2kikj — Akfkj + 6k 2 k 2 — Akikj + 2aklj + bujjlj + bujjk 
2ujiUjj — 2kikj — 4kfkj — 6k 2 k 2 — Akikj + 2al{lj + builj + bujjk 

(58) 

where i,j = 1, 2, 3, i < j. From the coefficient of e 3 , we obtain 

P(d){h} = -[A(l, 2)P(p 3 -p 2 - Pl ) + A(l, 3)P(p 2 - Pl - p 3 ) 

+ A(2, 3)P(p 1 - p 2 - p 3 )}e e ^ +e \ (59) 

Hence f 3 is in the form /3 = Be 6l+e2+03 where B is found as 

R 2)P(p 3 - Pl - p 2 ) + ,4(1, 3)P( P2 - Pl - gg) + A(2, 3)P(pi - p 2 - ga) 
±> = t . 

P{Pl +P2+ P3) 

(60) 

The coefficient of e 4 gives us 
e 2e 1+ e 2+ e 3 [jBp( ^ + ^ + A{1 ^ 2)A ^ 3)p{p2 _ ^ 

+ e ^ 29 ^ [BP( Pl + p 3 ) + A(l, 2)A(2, 3)P( Pl - p 3 )] 
+ e^ +e * +2e * [BP{ Pl + p 2 ) + 3)A(2, 3)P(pi - p 2 )] = (61) 

which is satisfied when 

£ = A(l, 2)A(1, 3)A(2, 3). (62) 
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The two expressions for B should be equivalent. This yields the three-Hirota 
solution condition (3HC) 



^2 P(<*iPi + a 2P2 + cr 3 p 3 )P(a 1 p 1 - cr 2 p 2 )P(a 2 p 2 - °?,P?)P{(?iPi ~ V3P3) = 0, 

(63) 



OY=±l 



r = 1, 2, 3 which turns out to be the below equation for eBo, 

(4a — b 2 )klklkl 2k\w 2 w 3 l 2 l 3 + 2k^wiw 3 lil 3 + 2k\wiw 2 lil 2 

+ 2k 2 k 3 w 2 w 3 l\ + 2kik 3 w\w 3 l 2 + 2kik 2 wiw 2 tl + 2k 2 k 3 w\l 2 l 3 

+ 2k\k 3 vo 2 H\l 3 + 2k\k 2 w\l\l 2 — 2k\k 3 w 2 w 3 l\l 2 — 2k\k 3 W\W 2 l 2 l 3 (64) 

- 2k\k 2 w 2 w 3 l\l 3 - 2k 1 k 2 w\w 3 l 2 l 3 - 2k 2 k 3 w\w 3 l\l 2 - 2k 2 k 3 w\w 2 l\l 3 

2 ,2,2 ;2 ,.2,2 ;2 ,.2,2 7.2„..2 ; 2 i,2„..2;2 ;2 ,.2,2 



12 Z/Z 7,Z Z/Z 7 Z Z/Z 7,Z Z/Z 7,Z Z/Z 7,Z Z/Z 

rv^t/Jg^ 1^2 3 2^1 3 2^3 1 3^2 1 3^1 2 



0. 



It is satisfied when a = 6 2 /4 or for some relations on k iy Wi and U which 
violate the solitonic property of the solution. The coefficients of e 5 and e 6 
vanish trivially. Similar to eKP, when a = 6 2 /4, eBo is integrable since it is 
transformable to Bo by the transformation 

u' — u, x' — x, y' = at + j3y, t' = t + py, 

where a = b/2 = y/a. When a 7^ o 2 /4, eBo is not integrable. The other 
relations which makes (3HC) satisfied are; 



Case 1. 


Any ki = 


0, i = 1,2,3, the rest are different, 


Case 2. 


ki — u}{, i 


= 1,2,3, 


Case 3. 


ki lif i 


= 1,2,3, 


Case 4. 


Ui = l h i 


= 1,2,3. 



By using any of these cases, we obtain the exact solutions of eBo. 
2.2.2 N = 4, Four-Hirota Solution of EBo 

Here we apply the Hirota direct method by using the ansatz which is used to 
construct four-soliton solutions. We take / = l+£fi+£ 2 f2+£ 3 j 3+£ 4 /4 where 
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f x = e e i+e e2 +e e3 + e 94 with 6 { = kiX+uJit+ky+ai for % = 1,2, 3, 4 and insert 
it into ([53)1 . We will only consider the coefficients of e m , m = 1,2,3,4,5 since 
the others vanish identically. By the coefficient of e 1 

P(D){l.ft + /x.l} = 2P{d){e ei + e 02 + e 03 + e di } = (65) 

we have the dispersion relation 

P( Vi ) = - kf - kj + alf + b^ih = (66) 

for i = 1, 2, 3, 4. From the coefficient of e 2 we obtain 

- P(d)f 2 = P(pi - p 2 )e 01+02 + P(pi - p 3 )e ei+03 + P(pi - p 4 )e 01+04 

+ P(P2 - Pz)e 02+03 + P( P 2 - pi)e 02+0A + P(p 3 - p 4 )e 03+04 . (67) 

Thus f 2 should be of the form 

f 2 = A(l, 2)e ei+d2 + A(l, 3)e 01+03 + A(l, A)e 01+0i 

+ A(2, 3)e 02+03 + A(2, A)e 02+0i + A(3, A)e 03+0 \ (68) 

where A(i,j) are obtained as 

A(i i) - P (P»-Pi) 
A{hj) ~ P(P, + P 3 ) 

2u)iUjj — 2kikj — Akfkj + Qkfk 2 — Akikj + 2aklj + bujilj + bujk 

2ujiL0j — 2k { kj — Akfkj — Qkfk 2 — Ak^ + 2al^j + buJilj + buOjU 

(69) 

for i, j = 1, 2, 3, 4 with i < j. The coefficient of e 3 gives 

(4) 

- P{d)h = Yl 3) p {v m -Pi- Pj) + A(i, m)P(p J - pi- p m ) 

i<j<m 

+ A{j, m)P(pi - Pj - p m )}e di+dj+6m (70) 

where (4) indicates the summation of all possible combinations of the four 
elements with i < j < m. Hence f 3 is of the form 

f 3 = 5(1, 2, 3)e 01+02+03 + 5(1, 2, A)e 01+02+0i 

+ 5(1, 3, A)e 01+03+0i + 5(2, 3, A)e 02+03+0 \ (71) 
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B(i,j,m) are obtained as 

A(i,j)P(p m -Pi-Pj) + A(i, m)P(pj - Pi - Pm) + A(J, m)P(pi - pj - p Ti 



B(i,j, m) 



P(Pi + Pj + Pm) 

(72) 

for i, j, m — 1, 2, 3, 4 with i < j < m. From the coefficient of e A we have 

P(D){f 4 A + / 3 ./i + h-h + fi-h + 1-/4 = 2P(d)U + 2P(L»){/i./ 3 } + P(D){f 2 .f 2 } 

= 0. 

(73) 

After simplifications we see that we should have 

B(i, j, m) = A(i, j)A(i, m)A( J} m) (74) 

for i,j,m = 1,2,3,4 with i < j < m. To be consistent the equations (|72j) 
and (|71jl should be equivalent. This gives us the condition 

ay=±l 



for i, j, m, r = 1, 2, 3, 4 with % < j < m, which becomes 
(Aa-b 2 )k 2 k 2 k 



(75) 



WiWjkL 

-\- 2,kjk m WjW m li -\- 2,kik m w iW m l j -\- 2,kikjWiWjl m -\- 2,kjk m w^ljl m 

+ 2k { k miUjlil m -\- 2kikjiu m lilj 2<kik m WjW m lilj 1kik m WiWjljl m (^6) 

*2kikjW jWrnljliYi *2kikjW iW m l jlm *2<kjk m WiW m lilj t 2kjk m WiWjlil m 

where m — 1, 2, 3, 4, i < j < m. Some of the cases except a = b 2 /4 which 
make this condition holds are; 

Case 1. Any two of ki = 0, % — 1, 2, 3, 4, the rest are different, 
Case 2. ki — u>i, % — 1, 2, 3, 4, 
Case 3. ki — l^i — 1, 2, 3, 4, 



20 



Case 4. Ui = l i} i = 1, 2, 3, 4. 

The equation remaining from the coefficient of e 4 is 

-P(d)U = e ei+92+e ^[B 123 P(p 4 -p 1 -p 2 -p 3 ) + B 12 ,P(p 3 ~p 1 -p 2 -p 4 )] 

+ B 1U P (p 2 -pi -p 3 -p 4 ) +#234 P (Pi -P3 -P4 ) + A(l , 2) A(3 , 4) P (pi +p 2 -p 3 -p 4 ) 

+ A(1, 3)A(2, 4)P(pi +p 3 -P2 -p 4 ) + 4)A(2, 3)P(p x +p 4 -p 2 -p 3 )] = 0. 

(77) 

Hence / 4 = C e ei+e2+e3+94 where C is obtained as 

C = -[A(l, 2)A(3, 4)P(pa + p 2 - p 3 - p 4 ) + 3)A(2, 4)P(p x + p 3 - P2 - p 4 ) 
+A(1, 4)A(2, 3)P(p a + p 4 - p 2 - Pa) + 5(1, 2, 3)P(p 4 - Pi - p 2 - p 3 ) 
+5(1, 2, 4)P(p 3 - Pi - Pa - Pa) + 3, 4)P(p 2 - Pl - p 3 - p 4 ) 

+P(2, 3, 4)P(p! -p 2 - p 3 - p 4 )] / P(pi + p 2 + Ps + p 4 )- 

(78) 

From the coefficient of e 5 we have 

2P(0)/ 4 + 2P(D){f 1 .f 3 } + P(P){/ 2 ./ 2 } = 0. (79) 
The simplifications give us that 

C = A(1, 2)A(1, 3)A(1, 4)A(2, 3)A(2, 4)A(3, 4). (80) 

To be consistent the equations (J78J) and (|50|l should be equal to each other. 
This yields the four-Hirota solution condition (4HC) 

A 

X)P(J><Pi) J] [Pfaft - oyp,)] = 0. (81) 
CTj=±l i=l 0<i<jr'<4 

In the hand, we know a case which satisfies both {3HC) and (AHC) auto- 
matically which is 

Case 1. Any two of hi = 0, i = 1, 2, 3, 4, the rest are different. 
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Here we give the graphs of the two- and four-Hirota solutions of eBo. We 
give arbitrary values to a, b, ki and w,i. From the dispersion relation, we 
obtain /j. We use these constants in the solutions and draw their graphs. 

i) N = 2, The Two-Hirota Solution of EKP: 



The constants are 

a = 5, b = 8, ki = 1, k 2 = 2, 




Figure 15: t=-2 Figure 16: t=2 
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Figure 17: t=4 Figure 18: t=6 

ii) iV = 4, The Four-Hirota Solution of EBo: 

The constants are chosen according to the Case 1 and the dispersion rela- 
tion. The constants are, 

a = 5, b = 8, ki = 0, k 2 = 0, k 3 = 2, k 4 = 3, 

w l = 5, w 2 = 4, w 3 = -6, w 4 = 1, 




Figure 19: t=-6 Figure 20: t=-4 
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Figure 21: t=-2 



Figure 22: t=2 




3 Conclusion 

In this work, we applied the Hirota direct method to non-integrable equa- 
tions. We have given two examples, the extended Kadomtsev-Petviashvili 
(eKP) and the extended Boussinesq equations. They are in general non- 
integrable equations. 

We have written bilinear and Hirota bilinear forms of these equations. 
Since the equations having Hirota bilinear forms automatically possess one- 
and two-Hirota solutions ( which have soliton-like behavior), we have focused 
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on three- and four-Hirota solutions. We have seen that both equations should 
satisfy a condition which we call three-Hirota solution condition (3HC) to 
have three-Hirota solution. While trying to obtain four-Hirota solutions of 
the equations we have come across another condition, four-Hirota solution 
condition (AHC). We have classes of solutions of these two conditions. 
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